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Abstract 


The  problem  of  geolocation  of  radar  emitters  has  been  studied  extensively  in  the  past 
decades.  Several  location  methods  such  as  the  methods  based  on  measuring  frequency 
difference  of  arrival  (FDOA),  time  difference  of  arrival  (TDOA),  angle  of  arrival 
(AOA),  etc.,  have  been  investigated  in  the  past.  However,  each  method  has  its  own 
advantages  and  limitations  in  terms  of  location  accuracy.  This  report  presents  a  multi¬ 
platform  geolocation  scheme  by  combining  TDOA  and  AOA  methods  for  three- 
dimensional  (3-D)  geolocation  of  radar  emitters.  A  mathematical  model  is  developed 
by  combining  sensor  information  such  as  AOA  measurements,  TDOA  measurements 
and  sensor  position  information  from  all  platforms.  A  least-squares  (LS)  solution 
based  on  the  model  is  derived  by  processing  the  sensor  information.  Through  the  use 
of  the  proposed  hybrid  TDOA/AOA  location  scheme,  location  ambiguity  due  to  using 
TDOA  method  alone  can  be  resolved  and  better  location  accuracy  can  be  achieved.  In 
addition,  unlike  single  platform  geolocation  techniques,  multi-platform  geolocation 
generally  results  in  quick  localization  of  radar  emitters.  Therefore,  it  will  reduce 
response  time  to  incoming  threats  and  improve  situation  awareness.  A  Monte-Carlo 
simulation  is  presented  using  additive  white  Gaussian  noise  to  demonstrate  the 
performance  of  the  proposed  hybrid  TDOA/AOA  geolocation  method. 


Resume 


Le  probleme  de  geolocation  d'emetteurs  radars  a  ete  etudie  de  fagon  extensive  au  cours 
des  decennies  precedentes.  Plusieurs  techniques  de  localisation,  telles  que  la  mesure  de 
la  difference  de  frequence  a  Tarrivee  (FDOA),  de  la  difference  de  temps  a  l’arrivee 
(TDOA)  et  de  Tangle  d’arrivee  (AOA),  ont  ete  investiguees  par  le  passe.  Cependant, 
chaque  methode  a  ses  avantages  propres,  ainsi  que  ses  limites  en  termes  de  precision 
en  ce  qui  regarde  la  localisation.  Ce  rapport  presente  une  approche  multi  plate-formes 
en  combinant  les  methodes  TDOA  et  AOA  pour  la  localisation  3-D  d'emetteurs  radar. 
Un  modele  mathematique  est  presente,  combinant  pour  toutes  les  plate-formes, 
l'information  provenant  des  detecteurs  telle  que  les  mesures  de  AOA,  TDOA  ainsi  que 
leur  position.  Une  solution  par  moindres  carres  (LS)  basee  sur  le  modele  est  derivee  en 
traitant  l'information  en  provenance  des  detecteurs.  En  utilisant  la  combinaison 
hybride  TDOA/AOA  proposee,  l’ambiguite  de  localisation  causee  par  TDOA 
seulement  peut  etre  resolue  et  une  plus  grande  precision  dans  la  localisation  spatiale 
peut  etre  atteinte.  Contrairement  a  la  localisation  basee  sur  une  seule  plate-forme,  la 
localisation  multi  plate-formes  fournira  des  mesures  rapides  et  precises  et,  par 
consequent,  ameliorera  la  connaissance  du  terrain  ainsi  que  le  temps  de  reponse  aux 
menaces.  Une  simulation  de  Monte-Carlo  est  presentee  pour  demontrer  la  performance 
de  la  technique  de  localisation  hybride  proposee.  L’ajout  de  bruit  blanc  gaussien  est 
utilise  durant  cette  simulation. 
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Executive  summary 


The  requirement  for  passive  geolocation  of  radar  emitters  remains  an  important 
problem  in  electronic  support  (ES)  for  maritime  and  littoral  surveillance.  Several 
location  techniques  such  as  the  ones  based  on  measuring  frequency  difference  of 
arrival  (FDOA),  time  difference  of  arrival  (TDOA),  angle  of  arrival  (AOA),  etc.,  have 
been  studied  extensively  in  the  past.  However,  in  general,  there  is  no  single  method 
that  will  provide  accurate  location  estimation  under  all  circumstances.  Each  method 
has  its  own  advantages  and  limitations  in  terms  of  location  accuracy. 

Both  TDOA  and  FDOA  location  methods  generally  provide  better  location  accuracy 
than  the  AOA-based  location  method  [2].  However,  the  TDOA-based  method  requires 
relatively  more  sensors  than  the  AOA-based  method  to  perform  geolocation  of  radar 
emitters.  For  instance,  the  TDOA-based  method  requires  at  least  three  properly 
distributed  sensors  for  two-dimensional  (2-D)  geolocation.  In  addition,  it  requires  that 
the  TDOA  measurements  at  all  sensors  be  performed  on  the  same  pulse  from  the  same 
pulse  train.  Otherwise,  location  ambiguity  will  be  introduced  if  the  TDOA 
measurements  are  taken  on  different  transmitted  pulses.  The  FDOA-based  method 
requires  relative  movement  between  the  transmitter  and  the  receiver.  The  AOA 
method,  on  the  other  hand,  requires  two  sensors  for  2-D  localization.  However,  it  is 
highly  range  dependent  and  a  small  error  in  the  angle  measurements  will  result  in  a 
large  location  error  if  the  sensor  is  far  away  from  the  transmitter.  Therefore,  to  achieve 
better  location  accuracy,  a  combination  of  two  or  more  location  schemes  should  be 
considered  in  order  to  complement  each  other. 

This  report  presents  a  multi-platform  geolocation  scheme  by  combining  TDOA  and 
AOA  location  methods  for  three-dimensional  (3-D)  geolocation  of  radar  emitters.  The 
proposed  hybrid  TDOA/AOA  location  scheme  is  achieved  by  utilizing  spatially 
separated  sensors  on  different  platforms,  such  as  UAVs  and  ship/land-based  platforms, 
to  improve  detection  range  and  location  accuracy.  The  radar  emitter  geolocation  is 
performed  by  processing  the  information  available  from  each  sensor,  including  AOA 
measurements,  TDOA  measurements  and  sensor  position  information.  Combining  the 
AOA  method  with  the  TDOA  method  can  assist  in  eliminating  location  ambiguity 
associated  with  TDOA  alone  and  can  enhance  location  accuracy.  In  addition,  the 
proposed  location  scheme,  unlike  single  platform  utilizing  individual  location  method, 
will  result  in  quick  and  better  location  estimates  of  radar  emitters.  Therefore,  it  will 
improve  ES  operation  and  situation  awareness. 
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The  theoretical  performance  and  effectiveness  of  the  proposed  location  scheme  is 
demonstrated  using  Monte-Carlo  simulation  with  additive  white  Gaussian  noise. 
However,  this  study  did  not  take  into  consideration  of  some  of  the  real-life  problems 
caused  by  multipath,  time  of  arrival  ambiguity  between  sensors,  asynchronous 
illumination  of  the  sensors  by  an  emitter  and  UAV  trajectory  selection.  Some  of  these 
issues  will  be  addressed  in  the  future. 


Du,  H.J.  and  Lee,  J.  P.Y.  2004.  Simulation  of  Multi-Platform  Geolocation  using  a  Hybrid 
TDOA/AOA  Method.  DRDC  Ottawa  TM  2004-256.  Defence  R&D  Canada  -  Ottawa. 


DRDC  Ottawa  TM  2004-256 


Sommaire 


L’exigence  de  la  geolocalisation  passive  d’emetteurs  radar  demeure  un  probleme 
important  dans  le  domaine  du  soutien  electronique  (SE)  de  la  surveillance  maritime  et 
cohere.  Par  le  passe,  plusieurs  techniques  de  localisation,  telles  que  celles  basees  sur  la 
mesure  de  la  difference  de  frequence  a  Tarrivee  (FDOA),  de  la  difference  de  temps  a 
Tarrivee  (TDOA)  et  de  Tangle  d’arrivee  (AOA),  ont  ete  etudiees  de  fagon  exhaustive. 
Toutefois,  en  regie  generate,  il  n’y  a  pas  de  methode  unique  qui  donne  une  estimation 
precise  de  la  position.  Chaque  methode  comporte  ses  propres  avantages  et  limitations 
en  termes  de  precision  de  la  localisation  en  toutes  circonstances.  Par  exemple,  la 
methode  TDOA  exige  en  general  au  moins  trois  capteurs  repartis  de  fagon  appropriee 
pour  la  localisation  bidimensionnelle  (2D)  avec  ambiguite  possible  de  la  localisation; 
la  methode  FDOA  exige  le  mouvement  relatif  entre  un  emetteur  et  un  recepteur.  Fes 
methodes  TDOA  et  FDOA  ont  generalement  une  meilleure  precision  que  le  systeme 
AOA.  Par  contre,  le  systeme  AOA  n’exige  que  deux  capteurs  pour  l’estimation  de  la 
position  en  2D.  Cependant,  il  est  extremement  sensible  a  la  distance,  et  une  faible 
erreur  des  mesures  d’ angle  donne  lieu  a  une  grande  erreur  de  localisation  si  le  capteur 
se  trouve  a  une  grande  distance  d’un  emetteur.  Pour  augmenter  la  precision  de  la 
localisation,  on  devrait  envisager  une  combinaison  de  deux  ou  plusieurs  systemes  de 
localisation  pour  obtenir  Tavantage  de  differentes  techniques  de  localisation 
complementaires. 

Dans  le  present  rapport,  nous  proposons  un  systeme  de  localisation  au  moyen  de  la 
combinaison  des  methodes  TDOA  et  AOA  pour  la  geolocalisation  tridimensionnelle 
(3D)  d’emetteurs  radar.  Fe  systeme  TDOA/ AOA  hybride  propose  est  realise  grace  a 
Tutilisation  de  capteurs  separes  dans  Tespace  (long  axe  de  reference)  sur  differentes 
plates-formes  pour  augmenter  la  portee  de  detection/couverture  et  ameliorer  la 
precision  des  mesures.  Nous  presentons  un  modele  mathematique  en  combinant 
Tinformation  des  capteurs,  telle  que  les  mesures  AOA,  les  mesures  TDOA  et 
Tinformation  de  position  des  capteurs  de  toutes  les  plates-formes.  Pour  localiser  les 
emetteurs  radar,  une  solution  des  moindres  carres  (FS)  basee  sur  le  modele  est  derivee 
par  le  traitement  de  Tinformation  des  capteurs.  Grace  a  la  combinaison  des  methodes 
TDOA  et  AOA,  on  peut  eliminer  T  ambiguite  de  la  localisation  due  a  la  TDOA  utilisee 
seule  et  obtenir  une  meilleure  precision  de  la  localisation  que  celle  de  TAOA  utilisee 
seule.  En  outre,  au  contraire  de  la  localisation  a  plate-forme  unique,  la  localisation  a 
plusieurs  plates-formes  produira  des  mesures  rapides  et  precises  et  ameliorera  done  la 
conscience  de  la  situation  et  le  temps  de  reponse  aux  menaces.  Nous  presentons  une 
simulation  Monte  Carlo  pour  faire  la  demonstration  de  la  technique  proposee  tout  en 
utilisant  du  bruit  blanc  gaussien  additif. 

Fa  performance  theorique  du  plan  propose  de  localisation  est  obtenue  au  moyen  d'une 
simulation  informatique  qui  prevoit  un  biais  et  un  bruit  additifs.  Il  est  continue  que  des 
estimations  plus  precises  de  femplacement  sont  obtenues  parce  que  le  plan  propose  de 
localisation  permet  une  estimation  et  le  retrait  subsequent  du  biais  des  mesures. 
Toutefois,  la  presente  etude  ifa  pas  tenu  compte  de  certains  des  problemes  reels  dus  a 
T ambiguite  causee  par  la  propagation  par  trajets  multiples  ou  les  ecarts  de  temps 
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d'arrivee  aux  capteurs,  a  rillumination  asynchrone  des  capteurs  par  un  emetteur  ou  a 
la  selection  de  la  trajectoire  des  VAT. 


Du,  H.  J.  and  Lee,  J.  P.  Y.  2004.  Simulation  of  Multi-Platform  Geolocation  using  a 
Hybrid  TDOA/AOA  Method.  DRDC  Ottawa  TM  2004-256.  R  &  D  pour  la  defense 
Canada  -  Ottawa. 
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IX 


1.  Introduction 


The  requirement  for  passive  geolocation  of  radar  emitters  remains  an  important 
problem  in  electronic  support  (ES)  for  maritime  and  littoral  surveillance.  Several 
location  techniques  have  been  extensively  studied  and  investigated  in  the  past 
including  angle  of  arrival  (AO A)  [1],  time  difference  of  arrival  (TDOA)  [2] [3]  and 
frequency  difference  of  arrival  (FDOA)  [4] [5] [6].  In  general,  there  is  no  single 
geolocation  method  that  will  give  accurate  location  estimation  under  all  circumstances. 
Each  location  technique  has  its  own  advantages  and  disadvantages.  The  TDOA 
method  is  based  on  measuring  the  time  difference  of  arrival  of  incoming  signals  at  two 
or  more  sensors.  The  TDOA  method,  in  general,  requires  at  least  three  properly 
distributed  sensors  for  two-dimensional  (2-D)  location  but  it  may  cause  location 
ambiguity  if  the  TDOA  measurements  are  performed  on  different  pulses  of  a  pulse 
train.  The  FDOA  method  is  based  on  measuring  the  frequency  difference  of  received 
signals  at  two  or  more  sensors  and  it  requires  relative  movement  between  sensors  and 
a  transmitter.  In  general,  both  TDOA  and  FDOA  have  better  location  accuracy  than 
the  AOA  method  [7].  The  AOA  method,  on  the  other  hand,  requires  only  two  sensors 
minimum  for  a  2-D  location  estimate.  However,  it  is  highly  range  dependent.  A  small 
error  in  the  angle  measurements  will  result  in  a  large  location  error  if  the  sensor  is  far 
away  from  a  transmitter.  To  achieve  better  location  accuracy,  a  combination  of  two  or 
more  location  methods  should  be  considered  in  order  to  complement  each  other. 

Fang  [8]  and  Schau  et  al  [9]  derived  a  closed-form  location  solution  when  the  number 
of  TDOA  measurements  is  equal  to  the  number  of  unknowns  (i.e.,  the  coordinates  of  a 
transmitter).  This  solution,  however,  cannot  make  use  of  extra  measurements  when 
they  become  available  to  improve  position  accuracy.  The  more  general  situation  with 
extra  measurements  by  utilizing  multiple  stationary  sensors  was  considered  by  Smith 
et  al  [10].  To  obtain  a  precise  position  estimate  with  reasonable  noise  levels,  the 
Taylor-series  method  [3][12]  is  commonly  employed.  Du  et  al  [13]  proposed  a 
Taylor-series  based  location  method  by  utilizing  one  or  two  moving  sensors.  The 
proposed  method  takes  into  account  both  random  noise  and  measurement  bias,  but  it 
requires  an  initial  guess  of  the  starting  point  to  start  iterative  location  estimation.  Chan 
et  al  [11]  proposed  an  explicit  approach  based  on  the  two-stage  least-squares 
computation  for  the  case  when  the  TDOA  estimation  error  is  small.  To  incorporate 
AO  As  for  a  better  position  estimate,  Cong  et  al  [7]  proposed  a  hybrid  TDOA/AOA 
scheme  for  2-D  mobile  user  localization  for  a  wireless  communication  application. 
Their  solution  is  based  on  a  two-stage  least-squares  algorithm  [11]  for  the  case  when 
the  TDOA  and  AOA  measurement  errors  are  small.  In  this  report,  we  propose  a  hybrid 
TDOA/AOA  location  scheme  which  extends  the  approach  in  [7]  for  a  2-D  location 
problem  to  3-D  geolocation  of  radar  emitters. 

The  proposed  hybrid  TDOA/AOA  scheme  is  achieved  through  combining  the  TDOA 
location  with  AOA  location.  As  the  TDOA  and  AOA  location  methods  complement 
each  other,  the  proposed  location  scheme  has  the  advantages  of  both  TDOA  and  AOA 
methods  to  achieve  better  location  accuracy.  The  proposed  location  scheme  utilizes 
spatially  separated  sensors  (long  baseline)  on  different  platforms  to  increase  detection 
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range/coverage  and  enhance  measurement  accuracy.  A  mathematical  model  is 
presented  by  combining  sensor  information,  such  as  AOAs,  TDOAs  and  sensor 
position  information  from  all  platforms.  A  least-squares  (LS)  location  solution  based 
on  the  model  is  derived  by  processing  the  information  available  from  each  sensor.  By 
the  combination  of  TDOA  and  AOA  methods,  the  location  ambiguity  associated  with 
the  TDOA  method  alone  can  be  resolved  and  better  location  accuracy  can  be  achieved. 
In  addition,  the  proposed  multi-platform  location  scheme,  unlike  single  platform  using 
individual  location  method,  will  result  in  quick  location  of  radar  emitters.  Therefore,  it 
will  reduce  response  time  to  incoming  threats  and  improve  situation  awareness.  The 
theoretical  performance  of  the  proposed  location  scheme  is  demonstrated  using  Monte- 
Carlo  simulation  with  additive  white  Gaussian  noise. 
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2.  Location  Methods  -  TDOA,  AOA  and  Hybrid 
TDOA/AOA 


2.1  The  TDOA  Method 

Fig.l  illustrates  the  geometry  of  TDOA  technique  using  multiple  stationary  sensors. 
The  time  of  arrival  (TO A)  is  defined  as  the  time  interval  for  signals  propagating  from 
the  transmitter  to  a  sensor,  as  shown  in  Fig.2.  The  TOA  is  dependent  on  transmitter¬ 
sensor  geometry  and  medium  characteristics.  In  the  case  of  a  constant  velocity 
medium  as  assumed  in  this  study,  the  time  interval  is  a  function  of  the  distance/range 
between  the  transmitter  and  the  sensor  [10].  Let  tt  be  the  time  interval  for  the  signal 
radiating  from  the  emitter  to  arrive  at  sensor  i  and  rt  be  the  distance/range  between  the 
emitter  and  the  sensor  /,  i.e., 


n  =  V<X  -  a)2  +  (ye  -  yt)2  +  <X  -  a)2 


(i) 


The  TDOA  between  sensors  i  and  j  is  then  defined  by  AtUj  as 


(2) 


where  c  is  the  speed  of  light,  ( xhyhz z)  is  the  position  coordinate  of  sensor  i  and  (xe,  ye, 
ze)  is  the  position  coordinate  of  the  emitter.  For  simplicity,  we  assume  that  the  TDOAs 
are  obtained  with  respect  to  the  first  sensor,  i.e.,  the  index  j  is  presumed  to  be  1  and  the 
index  i  is  presumed  to  run  from  2  to  k.  In  Eq.2,  the  Atitj  is  considered  as  the  error- free 
values  of  TDOAs.  In  practice,  however,  noise  is  always  present  in  the  TDOA 
measurements.  Therefore,  Eq.2  can  be  rewritten  as 


T.  —  T 

Att  j  =  - - 1  +  n 


'*i,i  ’ 


i  =  2,3,...,k, 


(3) 


where  n ,,  is  the  measurement  noise  associated  with  Att , .  In  order  to  solve  three 

unknowns  (xe,  ye,  ze),  a  minimum  of  three  TDOA  measurements  from  four  stationary 
sensors  is  required. 
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Sensor/:  (x,,  z) 


Figure  1.  The  geometry  of  TDOA  technique  using  multiple  stationary  sensors 
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► 
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Figure  2.  Time  of  arrivals  (TOA)  and  time  differences  of  arrival  (TDOA) 
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2.2  The  AOA  Method 


Fig.  3  shows  the  geometry  of  the  AOA  location  technique  in  the  Cartesian  coordinate. 
The  AOA  method  is  based  on  measuring  the  angles  of  arrivals  of  the  received  signals 
in  terms  of  azimuth  angle  9  and  elevation  angle  (j)  from  the  sensor  to  the  emitter.  The 
nonlinear  relationship  between  the  AO  As  and  the  sensor/emitter  coordinates  are 
expressed  by 


0{)  =  tan 


-i 


ye-ys 


(/)0  =  tan 


-i 


(4) 


where  0O  and  are  defined  as  the  error-free  values  of  the  azimuth  and  elevation 
angles,  ( xs,  ys,  zs)  are  the  sensor  position  coordinates,  and  (xe,  ye,  ze)  are  the  emitter 
location  coordinates.  In  practice,  however,  noise  is  always  present  in  the  AOA 
measurements.  Therefore,  Eq.4  can  be  rewritten  as 


6  =  80+  ne, 
</>  =  <f>o+n^ 


or 


ne  =0-0o. 


(5) 


(6) 


where  6  and  (j)  are  the  azimuth  and  elevation  measurements,  and  ne  and  n  are  the 
noise  related  to  the  azimuth  and  elevation  angle  measurements.  Taking  the  sine  of  both 
sides  of  Eq.6  and  multiplying  by  //  and  l2  respectively,  we  have 
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z 


Figure  3.  The  geometry  of  AO  A  location  technique  in  terms  of  the  azimuth  angle  6  and  the  elevation  angle  $  of 
received  signals  at  the  sensor 


lx  sin  n0  =  lx  sin(#  -  0O ), 
l2  sin  n ^  =  l2  sin(^  -  ), 


(7) 


where  lx  =  ^(xe  ~xs)2  +(ye  ~  ys)2  is  defined  as  the  length  of  x-j;  plane  projection  of 

vector  5E  and  l2  =  ^(xe  -x^)2  +(.ye  -  j^)2  +(ze  ~zs)2  is  the  distance  between  the 

emitter  and  the  sensor  in  the  Cartesian  coordinate  given  in  Fig.  3.  Considering  Fig. 3, 
we  have 
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(8) 


xe  -xs  =  /,  cos<90, 

ye  ~ys  =  A  sin^0» 

xe-xs  =  /2  cos  cos  6( 
ye-ys=l 2  COS  Sin  ^0>  (9) 

z,  -z,  =4  sin^o- 


Expanding  Eq.7  and  then  substituting  Eqs.  8  and  9  into  Eq.7  gives 


/jSin^  =(x  -x  )sin0-(>>  )cos0, 


l2s\nn4=(xe-xs) 


sin^ 

cos<9„ 


-(ze-zs)cos</>. 


(10) 


To  simplify  the  formulation  of  the  problem,  we  assume  that  AO  A  measurement  noise 
are  very  small,  i.e.,  |  ne  |  and  |  n  \  «  1  and  therefore,  sin  (ne)  ~  ne,  sin  (n  )  ~  n  and 
cos  (00)  ~  cos  (9).  Based  on  these  assumptions,  we  can  approximately  rewrite  Eq.10  in 
a  linear  form  as 


0  S  ~(xe  -  xjsin#  +  (ye  - ys) cos#  +  /,  »*, 

0«-(xe-xJ-^^  +  (ze-zs)cos^  +  /2  nr 
cos  0 


or 


q-.  0>e->Qcosfl 


A 


A 


0  -.  (^-^)sin^  (  (ze-z,)cos^  |  ^ 

l2COS0  l2  * 


(11) 


(12) 


As  can  be  seen  from  Fig. 3,  Eq.12  is  derived  based  on  the  angle  measurements  from  a 
single  sensor.  In  the  case  of  multiple  sensors,  e.g.,  k  sensors,  Eq.12  can  be  rewritten  as 
a  set  of  equations  as 
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|  (ye-y„j) cos  <9, 

^ n,\  n,\ 

(Xe-Xn,JSin^  ;  (Ze  -Z«,,)C0S^ 

4,2  cos  0  /„,2 


n  =  \,...,k, 


where  the  subscript  n  in  Eq.  13  denotes  the  n  sensor. 


(13) 


2.3  The  Hybrid  TDOA/AOA  Method 

Combining  Eq.3  and  Eq.13,  and  then  re-arranging  it  in  a  matrix  form  with  i  =  2 
and  n=  1 ,2 . k.  we  have 


m  =  g(X )  +  n  , 


(14) 


where  X  =  [xe  r  ze  ]r  is  the  vector  of  the  unknown  variables  and 


^2,1 

^3,1 

1 

sc- 

1 

i _ 

(r3  _  ri)/ C 

••  ^  ©  O 

(rk-r\)/c 

,g(X)  = 

~  (xe  - xi,s) sin 6*1  //|,|  +  (ye  -  J;i,Jcos6»,  //u 
-  (xe  -xi,,) sin  </>i  /(A.2  cos 6> )  +  (ze  -  zM )  < cos  ^  / . /u 

_  0 

-  (*«  -  ) sin  <4  /  4,i  +  (34  -  34, s )  cos  ek  /  /w 

(xe  -  xi,s)sin  A-  /(4,2  cos^)  +  (ze  -  zks)co%</>k  !lkl_ 
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and  n  is  the  vector  of  TDOA  and  AOA  measurement  noises,  i.e., 


n=[n 


'Atj 


Atk 


net  n 4 


(16) 


In  the  following,  an  LS  solution  based  on  the  combination  of  TDOA  and  AOA 
measurements  is  introduced  to  solve  this  nonlinear  problem. 
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3.  Multi-Platform  Geolocation  Using  the  Hybrid 
TDOA/AOA  Method 


3.1  Geolocation  from  UAV  and  Ship/Land-Based  Sensors 


Unmanned  Aerial  Vehicles  (UAVs)  today  are  gaining  wide  acceptance  as  valuable 
reconnaissance,  surveillance,  targeting  and  intelligence  gathering  tools  to  support 
various  missions.  However,  one  of  the  common  uses  for  a  UAV  today  is  to  find  a 
target  and  determine  its  location.  Fig.  4  illustrates  a  possible  operational  scenario  of 
passive  geolocation  from  UAV  and  ship/land-based  sensors  using  the  proposed  hybrid 
TDOA/AOA  method. 

3.2  Formulation  of  the  Problem 

Instead  of  using  multiple  stationary  sensors,  Fig.4  shows  the  geolocation  scenario 
where  two  time-synchronized  ESM  sensors  are  utilized  with  sensor  1  on  the  ship/land- 
based  platform  and  sensor  2  on  the  UAV  platform.  The  multiple  TDOA  measurements 
are  collected  from  the  two  sensors  over  a  time  period  along  the  UAV  path.  The 
multiple  AOA  measurements  (i.e.,  elevation  and  azimuth  angles)  are  obtained  from  the 
ship/land-based  sensor. 

For  simplicity  of  analyzing  the  problem,  Fig.4  is  re-illustrated  in  Fig. 5  with  sensor  1 
being  the  origin  (i.e.,  x  =  0,  y  =  0  and  z  =  0)  of  the  coordinate.  The  coordinates  of 
sensor  2  are  defined  by  (x^jy 2^,2)  at  the  zth  time  instant  along  the  UAV  trajectory. 

•  Multiple  TDOA  Equations 

Eq.14  defines  a  set  of  nonlinear  equations  as  a  function  of  unknown  variables  (xe,ye, 
ze).  Since  solving  these  nonlinear  equations  is  difficult,  an  alternative  [11]  is  to 
transform  these  equations  into  another  set  of  linear  equations  through  introducing 
another  unknown  variable,  rx. 

For  simplicity  of  the  following  derivation,  re-writing  Eq.3  without  consideration  of 
TDOA  measurement  noise  gives 


c 


(17) 


i  =  2,3,...,k, 


or 


10 


DRDC  Ottawa  TM  2004-256 


Sensor  1 

Ship/Land-Based 

Platform 


Sensor  2 


Mini-UAV 

tdoa=toa,-toa2 


Figure  4.  Geolocation  using  the  hybrid  TOO  A/ AO  A  method  from  a  mini-UAV  in  conjunction  with  a 
ship/land-based  platform 


Figure  5.  Geometry  between  an  emitter  and  two  sensors. 
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(18) 


cAti,i  =ri~rl, 

i  =  2,3,...,k, 

where  rx  is  defined  as  the  range  between  sensor  1  and  the  emitter,  and  rt  is  the  range 
between  sensor  2  and  the  emitter  at  the  /th  time  instant  where  i  ^  1  (i  =  2,...,  k)  to 
distinguish  rt  from  rx. 

In  the  following,  we  define  ri>x  as  the  range  difference  between  rt  and  rx  at  the  ith  time 
instant  along  the  UAV  path.  Therefore,  Eq.18  can  be  rewritten  as 


rui=cAtui=ri-ri’ 

i  =  2,3 


(19) 


Re-arranging  Eq.19  and  then  taking  the  square  of  both  sides  of  the  equation  gives 


r-  ={rix+rxf, 

i  =  2,3  ,...,k. 


(20) 


As  can  be  seen  from  Fig.5,  we  have 


and 


2  ,  2  ,  2 

=  +ye  +Ze 


(21) 


r,2  =  (xi,  2  -  xef  +  (yU2  -  yef  +  (Zi,  2  -  Ze)2’ 
i  =  2,3  ,...,k. 


(22) 
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where  (.y,2-  v,.2  zu2)  are  the  coordinates  of  sensor  2  at  the  time  instant  along  the  UAV 
trajectory.  Substituting  Eqs.21  and  22  into  Eq.20  gives 


2 


C n/  -Ki)  =  ~xi,2xe  -yu 2ye 

i  =  2,3,...k, 


■Zi,2Ze 


ri,A’ 


(23) 


where 


Ki  ~  Xi, 2  +  »,2  +  Zi, 2  • 


(24) 


Rewriting  Eq.23  with  i  =  2.3 . k  in  a  matrix  form  gives 


=  gjX  , 


(25) 


where 


<N 

1 _ 

-k2' 

X2,2 

yi,2 

Z2,2 

1 - 

1 

3 

Ki 

II 

K>|  - 

2 

r3,l 

~K, 

»  gi  = 

~  X3,2 

-y3, 2 

"Z3,2 

-r34 

2 

Lr*.i 

l 

1 

<N 

1 

_ 1 

-yk,  i 

~Zh,2 

- 1 

1 

(26) 


and  X  =  [xe  rx  ]T  is  the  vector  of  unknown  variables  xe,  ye,  ze  and  ri. 

As  can  be  seen,  Eq.  25  is  derived  without  consideration  of  TDOA  noise.  In  practice, 
however,  noise  is  always  present  in  the  measurements.  To  have  better  formulation  of 
the  problem,  an  error  vector  (p 7  is  introduced  into  Eq.25  as 


mi = gix + <pi  ■ 


(27) 


DRDC  Ottawa  TM  2004-256 


13 


Multiple  AOA  Equations 


Considering  Eq.13  with  sensor  1  at  the  coordinate  origin  as  shown  in  Fig. 5,  i.e.,  (xi>s, 
yUs,  zUs)=  (0,  0,  0),  we  have  AOA  equations  as 


0  =  -xe  sin  0i  +  ye  cos  0i  +  li  X  ne  , 

0  =  -Xe  SUU^  +  Ze  C0S  4>i  +  li,2  \  , 

cos  0. 


(28) 


where  0t  and  (/).  are  the  azimuth  and  elevation  angles  measured  at  sensor  1  at  the  zlh 
time  instant.  Eq.28  can  be  rewritten  in  a  matrix  form  with  i  =  2,3 ,...,£  as 


m2=g2X  +  <P2’ 


(29) 


where 

*=ix  ^  rJT > 


(30) 


is  the  vector  of  unknown  variables  and 


^2  —  [4,1%,  h,2n(/>2  ’  *  *  h,inek  h,2nfa  ]  ’ 


(31) 


is  the  error  vector  associated  with  the  AOA  equations,  and 
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0 

0 

0 


-  sin  02 

-sin  (j)2/  cos  d2 


cos  e2  o  o 

0  cos  (f>2  o 


g 2  = 


-  sin  6k 


-sin  (j)J  cos  0k 


cos  6k 
0 


0 

cos  </>k 


0 

0 


(32) 


•  Multiple  Hybrid  TDOA/AOA  Equations 


Combining  Eq.25  and  Eq.29  gives 


ntj 

gi~ 

x  + 

<p,~ 

— 

m2 

g2_ 

_(P2_ 

(33) 


Eq.33  represents  a  set  of  overdetermined  location  equations  as  the  number  of 
measurements  is  greater  than  the  number  of  unknowns.  In  the  presence  of  noise,  this 
set  of  equations  will  not  meet  at  the  same  point  and  the  proper  solution  is  the  ( xe ,  ye,  ze, 
r\ )  that  best  fit  these  equations. 

3.3  The  Least-Squares  Approach 


The  vector  X  of  unknown  variables  xe,y(,  ze  and  r\  can  then  be  estimated  by  the 
following  weighted  LS  calculation 


X  =  argmin 

WU2 

ntj 

_ 

gi 

X 

V 

m  2 

_g  2  _ 

J 

f  r 

T 

r-  -T 

y1 

r  -| 

T 

r 

gi 

w 

gl 

gi 

w 

ntj 

w 

-§2~ 

/ 

gi_ 

m2 

(34) 


(35) 
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The  emitter  position  coordinates  xe,  ye  and  ze  can  be  picked  out  by  pre-multiplying  the 
both  sides  of  Eq.35  by  a  3x4  matrix  composed  mostly  of  zeroes  as 


~1 

0 

0 

o] 

T 

' 

y7 

_ 

T 

ye 

~ 

0 

1 

0 

0 

gl 

w -1 

gl 

gi 

W  1 

nij 

?e_ 

0 

0 

1 

oj 

_g 2  _ 

) 

_g  2  _ 

m2 

(36) 


Note  W in  Eq.36  is  defined  as  a  weighting  matrix.  In  the  case  where  q>\  and  qy  in  Eq.33 
are  Gaussian  random  vectors  and  W is  suitably  chosen  as  the  covariance  matrix  of  [q>\ 
T  qyT]T,  i.e.,  W=  E([^i T  qyT]T[q>xT  ^T]),  then  the  position  estimation  in  Eq.36  becomes 
a  maximum-likelihood  (ML)  estimate  [11].  However,  the  weighting  matrix  required  in 
LS  is  not  easy  to  determine  and  therefore  it  is  not  discussed  in  this  report.  To  simplify 
the  problem,  we  define  the  weighing  matrix  W=  I  (an  identity  matrix)  here. 

The  LS  approach  to  solve  the  nonlinear  problem  is  based  on  the  assumption  of 
independence  among  xe,ye,  ze  and  rx.  They  are  then  solved  by  the  LS  given  in  Eq.36. 
An  improved  estimate  of  the  unknown  variables  can  be  obtained  via  another  LS 
computation  [11]  by  applying  the  computed  result  (xe,ye,ze)  to  the  known 
relationship  in  Eq.21  to  solve  rx. 
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4.  Monte-Carlo  Simulation 


The  performance  of  the  proposed  location  scheme  has  been  evaluated  using  computer 
simulation.  One  type  of  computer  simulation  is  Monte-Carlo  simulation.  Monte-Carlo 
simulation  is  a  stochastic  technique  used  to  investigate  and  solve  mathematical 
problems  where  uncertain  or  random  variables  are  presented.  The  term  "Monte  Carlo” 
comes  from  the  name  of  a  city  in  Monaco  [14]  where  main  attractions  are  casinos 
running  games  by  exploiting  the  random  behaviour  of  each  game.  Monte-Carlo 
methods  were  originally  developed  for  the  Manhattan  Project  during  World  War  II 
because  of  the  similarity  of  statistical  simulation  to  games  of  chance.  However,  they 
are  now  applied  to  a  wide  range  of  problems. 

In  the  following,  Monte-Carlo  simulation  was  used  to  evaluate  the  performance  of  the 
proposed  geolocation  scheme  with  the  consideration  of  TDOA  and  AOA  measurement 
noise.  The  measurement  noise  for  TDOAs  and  AOAs  were  assumed  to  be  white  noise 
with  Gaussian  distribution.  During  each  run  of  the  Monte-Carlo  simulation,  random 
measurement  noise  was  generated  based  on  its  defined  statistical  distribution  (such  as 
normal  distribution,  mean  and  variance)  and  then  emitter  position  was  estimated  based 
on  this  randomly  generated  noise. 

In  the  following,  the  simulation  was  performed  based  on  the  geolocation  scenarios 
given  in  Figs.4  and  5.  For  each  simulation  scenario,  10,000  simulation  runs  were 
performed  and  the  emitter  location  estimates  (xe9ye9ze)  were  then  computed  by 
averaging  the  10,000  location  estimate  results  over  the  10,000  simulation  runs. 

To  illustrate  the  accuracy  of  location  estimation,  the  location  estimation  error  is  used 
and  defined  by 


Location  Estimation  Error  =  ^(xe  -xe)2  +  (ye  -ye )2  +  (ze  ~ze)2 


(37) 


where  (xe,ye,ze)  are  the  location  estimates  of  the  emitter  and  ( xe9ye9ze )  are  the 
emitter  true  positions. 
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Simulation  scenario:  geolocation  from  a  ship/land-based  sensor  accompanied  by  a 
mini-UAV  based  sensor 

In  the  following,  Figs.  6(a)  and  7(a)  show  the  geometry  between  the  emitter  and  the 
sensors  as  well  as  the  UAV  sensor  trajectory.  We  assume  that  the  emitter’s  true 
position  is  at  (xe,  ye,  ze)=(  10.0,  10.0,  2.0)  kilometers  (km)  in  the  coordinate  given  in 
Figs.  6(a)  and  7(a).  As  can  be  seen  in  Figs  6(a)  and  7(a),  the  crosses  represent  the 
positions  of  the  UAV  sensor  at  different  position  along  the  UAV  fight  trajectory  and 
the  square  at  the  origin  of  the  coordinate  represents  the  position  of  the  ship/land-based 
sensor. 

Based  on  the  given  geometries  in  Figs  6(a)  and  7(a),  the  distance  between  sensor  1  and 
the  emitter  is  about  14.3  km,  while  the  distance  from  sensor  2  to  the  emitter  varies 
from  3.3  km  to  4.2  km.  The  associated  range  difference  varies  from  10.1  km  to  11.0 
km,  corresponding  to  TDOA  variation  from  33.5  to  36.5  microseconds  (psec).  The 
azimuth  and  elevation  angles  from  sensor  1  to  the  emitter  are  about  0.785  radians  (i.e., 
about  45  degrees)  and  0.141  radians  (i.e.,  about  8.1  degrees)  respectively. 

To  evaluate  the  performance  of  the  proposed  method  in  handling  measurement  noise, 
the  noisy  measurements  were  produced  by  adding  white  Gaussian  noise  into  the 
TDOA  and  AO  A  measurements.  Figs.  6  and  7  show  location  estimation  based  on 
noisy  TDOA  and  AO  A  measurements.  The  simulated  standard  deviation  of  the  TDOA 
measurements  was  10.0  and  30.0  nanoseconds  (nsec)  respectively,  while  the  simulated 
standard  deviation  for  the  AOA  measurements  was  fixed  at  0.2  degrees.  The  location 
estimation  shown  in  Figs.  6  and  7  were  performed  by  collecting  six  TDOA  and  one 
AOA  measurements  for  each  estimation  interval.  The  associated  estimation  errors  at 
every  estimation  interval  are  given  in  Figs.  6(b)  and  7(b).  To  compare  the  simulation 
results  at  different  TDOA  noise  levels,  the  estimation  errors  given  in  Figs.  6(b)  and 
7(a)  were  re-plotted  and  given  in  Fig.  8.  As  can  be  seen,  an  increase  in  TDOA  noise 
level  increases  estimation  error,  as  we  would  expect. 

A  numerical  simulation  with  increased  number  of  TDOA  measurements  was  also 
performed  and  the  associated  estimation  results  are  given  in  Fig.  9.  As  can  be  seen,  an 
increase  in  the  number  of  TDOA  measurements  improves  estimation  accuracy. 
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(a) 


(b) 

Figure  6. .  Monte-Carlo  (10,000  runs)  location  estimation  with  TDOA’s  std  dev=10.0  nsec  and  AO  A’ s  std  dev=0.2 
deg.  (a)  Location  estimates  against  emitter  true  position;  (b)  Location  estimation  error. 
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(b) 


Figure  7.  Monte-Carlo  (10,000  runs)  location  estimation  with  TDOA’s  std  dev=30.0  nsec  and  AO  A’ s  std  dev=0.2 
deg.  (a)  Location  estimates  against  emitter  true  position;  (b)  Location  estimation  error. 
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Location  Estimation  Error  in  Distance  (km) 


Figure  8.  Monte-Carlo  (10,000  runs)  location  estimation  errors  with  TDOA’s  std  dev -10.0  and  30.0  nsec 
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Location  Estimation  Error  in  Distance  (km) 


Figure  9.  Monte-Carlo  (10,000  runs)  location  estimation  errors  with  different  number  of  TOO  A 
measurements 
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5.  Summary  and  Conclusions 


This  report  proposes  the  development  of  a  multi-platform  location  technique  using  a 
hybrid  TDOA/AOA  scheme  for  passive  radar  emitter  geolocation  applications. 
Implicit  in  the  development  are  the  assumptions  that  the  mechanisms  for  measuring 
TDOAs  and  AOAs  on  radar  pulses  exist,  and  that  TDOA  and  AOA  measurement 
errors  are  small.  By  incorporating  a  set  of  the  AOA  equations,  we  extend  the  previous 
TDOA  only  location  estimators  to  solve  the  hybrid  TDOA/AOA  equations  for  3-D 
location  applications.  By  introducing  an  intermediate  variable  ru  the  nonlinear 
equations  related  to  position  estimates  are  transformed  into  a  set  of  the  linear  equations 
which  are  a  function  of  the  unknown  variables  and  the  intermediate  variable.  Monte- 
Carlo  simulation  results  demonstrate  that  the  proposed  location  scheme  is  effective 
and  that  it  converges  fast  with  no  need  of  an  initial  guess  of  emitter  position.  However, 
further  studies  need  to  be  done  in  the  future  to  investigate  the  performance  of  the 
proposed  scheme  in  handling  higher  TDOA  and  AOA  noise.  Future  studies  also 
include  applying  another  LS  computation  by  applying  the  computed  result 
(xe, ye9ze9rl)  to  the  known  relationship  to  obtain  an  improved  position  estimate 

(xe,ye,ze)  as  well  as  incorporating  TDOA/AOA  measurement  variances  as  a 
weighting  matrix  in  the  LS  to  achieve  a  maximum-likelihood  (ML)  position  estimate. 
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